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20 ABSTRACT: A tropospheric ozone Differential Absorption Lidar (DIAL) system, developed 

21 jointly by the University of Alabama at Huntsville and NASA, is making regular observations of 

22 ozone vertical distributions between 1 and 8 km with two receivers under both daytime and 

23 nighttime conditions using lasers at 285 and 291 nm. This paper describes the lidar system and 

24 analysis technique with some measurement examples. An iterative aerosol correction procedure 

25 reduces the retrieval error arising from differential aerosol backscatter in the lower troposphere. 

26 Lidar observations with coincident ozonesonde flights demonstrate that the retrieval accuracy 

27 ranges from better than 10% below 4 km to better than 20% below 8 km with 750-m vertical 

28 resolution and 10-min temporal integration. 

29 1 . Introduction 

30 Ozone is a key trace-gas species within the troposphere. Although ozone is a precursor of 

31 the hydroxyl radical, which reacts with most trace species in the atmosphere, ozone is also a strong 

32 greenhouse gas influencing the climate by its radiative forcing and is a direct pollutant [1]. In situ 

33 photochemistry and dynamic processes largely govern the distribution of tropospheric ozone [2], 

34 Measuring ozone variability at high spatial and temporal resolution increases our understanding of 

35 tropospheric chemistry [3,4], planetary boundary layer (PBL)-free-tropospheric exchange [5, 6] 

36 stratosphere-troposphere exchange (STE) [7-9] and the impact of lightning-generated NOx on 

37 tropospheric ozone [10-13], 

38 Several techniques currently exist for making range-resolved measurements of 

39 tropospheric ozone. The most co mm on technique is the balloon-borne electrochemical 

40 concentration cell (ECC), which has monitored ozone since the 1960’s. The ozonesonde profiles 

41 ozone with a 100-m vertical resolution from the surface to 3 5 -km altitude with the accuracy of 
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42 5-10% [14, 15]. Ozonesondes are attractive because of their low upfront cost and 

43 well-characterized behavior. However, they are not suitable for making continuous measurements 

44 because of logistical considerations. Interesting atmospheric phenomena that vary over periods 

45 less than one day are particularly difficult to monitor using balloon ozonesondes. Satellite 

46 observations can derive total column ozone [16], stratospheric ozone [17-21], and extend 

47 measurements to altitudes that are inaccessible to ozonesondes. More recently, high-quality 

48 satellite observations of tropospheric ozone are becoming available [17, 22-30]. Although the 

49 satellite measurements can produce global maps of ozone, their current measurement uncertainties 

50 along with their coarse spatial and temporal resolution limit their ability to observe short-term 

51 variations in ozone. Lidars can supplement these techniques when a requirement exists for ozone 

52 retrievals with higher temporal (from 1 min to several hours) and vertical resolution (from tens of 

53 meters to 2 km). For example, lidars of the Network for the Detection of Atmospheric 

54 Composition Change (NDACC) [31, 32] are providing long-term observations of ozone as well as 

55 aerosol, temperature, and water vapor. Although the upfront costs are considerably higher than for 

56 a balloon ozonesonde operation, lidars can acquire profiles continuously under both daytime and 

57 nighttime conditions. The spatial and temporal resolution of a lidar is more than sufficient to 

58 characterize short-term ozone variations for the photochemical studies of vertical processes. 

59 Differential Absorption Lidar (DIAL) has been successfully used to measure ozone within the 

60 PBL [33, 34], the free-troposphere [35-40], and the stratosphere [41-44] for several decades. DIAL 

61 is evolving from ground-based and airborne systems to systems suitable for long-term deployment 

62 in space [45], The technique derives ozone concentrations by analyzing how rapidly the 

63 backscattered signals at two separate but closely spaced wavelengths, one strongly absorbed by 

64 ozone and the other less strongly absorbed, diminish with altitude. This measurement does not 
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65 require the absolute signal intensities but only how the relative change of the two signals with 

66 respect to altitude. The wavelengths optimize the differential absorption of ozone and minimize 

67 differential extinction due to aerosols, SO 2 , and other species. Using electronically-gated detection 

68 permits range-resolved measurements to a resolution as small as several meters over acquisition 

69 times of several minutes. The location of the ozone DIAL discussed in this paper in the 

70 southeastern United States provides a unique observational site within an interesting scientific 

71 area [46] to study trace gas transport at the mid-latitudes for both the polluted PBL and the free 

72 troposphere. 

73 2. System Description 

74 Housed in the Regional Atmospheric Profiling Center for Discovery (RAPCD), the 

75 tropospheric ozone DIAL system is located at 34.7250°N, 86.6450°W on the campus of the 

76 University of Alabama in Huntsville (UAHuntsville) within the Huntsville city limits at an 

77 elevation of 206 m ASL. It is designed for measurements within the PBL and the free-troposphere 

78 during both daytime and nighttime. Because of UAHuntsville ’s location and occasional high 

79 temperature and humidity conditions, heavy aerosol pollution is sometimes present. Compared 

80 with the clean free-troposphere, these aerosols require a larger dynamic range for the detection 

81 system because of their larger optical depth. Moreover, the rapid change of aerosol concentrations 

82 (e.g., due to convective activity) increases the measurement uncertainty for DIAL within the PBL 

83 and lower troposphere. Judicious system-design choices and an effective aerosol-correction 

84 scheme allow this system to produce high-quality ozone profiles under a variety of conditions. 

85 2. 1 Wavelength selection 

86 The selection of the 285- and 291-nm wavelengths results from the balance of the 
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87 following three considerations: 1) optimizing the altitude range to make retrievals; 2) reducing the 

88 impact of the solar background during daytime operation; 3) reducing the impact of aerosol 

89 interference upon the ozone retrieval. The DIAL wavelength selection is flexible and optimized 

90 for the local ozone distribution, the absorption arising from non-ozone species, the measurement 

91 range and the specific system configuration including the output power, telescope mirror size, and 

92 the photomultiplier’s (PMT’s) dynamic range. Numerous publications (e.g., Megie et al. [47]) 

93 discussed the optimum wavelengths for tropospheric systems. Although shorter wavelengths can 

94 provide higher measurement sensitivity arising from the larger ozone differential cross-section, 

95 they limit the maximum measurable range due to stronger attenuation of ozone absorption and 

96 Rayleigh (molecular) extinction, and thus, require more signal acquisition time. In addition, the 

97 shorter wavelengths require more dynamic range of the detection system and might require more 

98 altitude channels. With the current transmitter power, the on-line wavelength of 285 nm allows us 

99 to measure ozone up to 9 km under a clear sky and 7 km under aerosol loading with a 10-min 

1 00 temporal resolution. Because of the increasing solar background during daytime operations, we 

101 choose 291 nm as the off-line wavelength. Longer wavelengths will cause a significant increase in 

1 02 the solar background and reduce the signal-to-background ratio. To measure both wavelength 

1 03 channels using the same PMT and simplify the system design, we used a band-pass filter with a 

1 04 central wavelength of 286.4 nm and a full width at half maximum (FWHM) of 1 1 nm whose 

105 transmittance is < 10" 8 at wavelengths longer than 300 nm. For a band-pass filter, the integrated 

1 06 sky background over the filter bandwidth and the dark counts actually determine the background 

107 for both off-line and on-line wavelengths. For our lidar configuration, the 285 and 291-nm 

1 08 wavelength region can provide sufficient signal-to-background ratios at 8 km under most sky 

1 09 conditions. The retrieval errors due to aerosol interference are a concern in the PBL and lower 
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110 troposphere. These errors are not a simple function of the wavelength separations because 

111 reducing the separation to reduce the aerosol differential backscattering will also decrease the 

112 differential ozone cross-section. These errors are sensitive to the local aerosol composition, size 

113 distribution, and vertical profile. Although the aerosol interference can be lower when our on-line 

114 wavelength extends to the steepest part of the ozone absorption cross section, this will 

115 significantly sacrifice the maximum measurable range. Therefore, the 285-291 nm pair is the 

116 optimal choice to balance the maximum measurable altitude, the impact of aerosol differential 

117 backscattering, and the impact of solar background. 

118 2.2 Hardware components 

119 Table 1 lists the characteristics of the RAPCD ozone DIAL system. The transmitter 

120 consists of two identical dye lasers pumped by two separate, frequency-doubled Nd:YAG lasers 

121 (Figure 1). A pulse generator triggers each laser pulse with a 25-ms separation between the 

1 22 alternate pulses. The dye lasers are software-controlled to select the user-defined wavelength. The 

1 23 knife-edge method [48] determines that the divergences of both UV laser beams are less than 1 

124 mrad. A 0.75-m triple grating monochromator (Acton Research Corporation) indicates that the 

1 25 actual wavelengths of the outgoing UV lasers are 285 and 291 nm within an uncertainty of 0. 1 nm. 

1 26 The receiving system currently operates with two separate telescopes as shown in Figure 

1 27 2. The high-altitude receiver uses a 40-cm Newtonian telescope, and the low-altitude channel 

128 employs a 1 0-cm Cassegrain telescope. The large telescope system routinely makes measurements 

129 from 3 to 8 km and on occasion has reached 12 km. Employing a 1.5-mrad field-of-view (FOV), 

1 30 the large telescope achieves full overlap between the laser and receiver at about 3 km. Larger 

131 FOVs lower the altitude at which full overlap occur but significantly increase solar background. 

1 32 The small telescope system currently retrieves ozone between 1 and about 5 km with a typical 
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1 33 FOV of 4.3 mrad. The future plan is to extend the retrievals down to about 200 m with an 

1 34 additional altitude channel in the small telescope. The band-pass filters used to restrict the solar 

1 35 background for both receivers have a transmittance of 35% at 285 nm and 20% at 291 nm. 

1 36 The detection system of the RAPCD ozone DIAL uses both photon counting and analog 

1 37 detection to facilitate operations over both altitude channels. This detection combination provides 

138 the linearity of the analog signal in the strong-signal region and high sensitivity of the 

139 photon-counting signal in the weak-signal region. A EMI 9813 QA PMT, which has been used 

140 extensively for many years on a number of Goddard Space Flight Center lidar system [49, 50] is 

141 used in the high-altitude channel while a small Hamamatsu 7400 PMT is used in the low-altitude 

142 channel. A photodiode detects the outgoing laser pulses, which trigger both the PMT gating 

143 circuits and the LICEL Transient Recorder (TR40-80, LICEL Company, Germany). The 

144 LICEL-TR offers the advantage of increased dynamic range by providing simultaneous 

145 measurements using both analog detection and photon counting. The LICEL-TR’ s highest 

146 temporal resolution is 25 ns corresponding to a fundamental range resolution of 3.75 m. It is 

147 necessary to gate the high-altitude channel off for the first 10-15 ps and the low-altitude channel 

148 for the first 1 ps to maintain the PMT’s linearity and minimize the impact of signal-induced bias 

1 49 on the background count rate. 

1 50 3. Data processing 

151 3.1 Raw data processing 

152 Several operations, designed to improve the measurement precision, occur before the 

1 53 ozone retrieval. First, average the signal returns over 10 min and 150 m. The temporal resolution 

1 54 of the retrieval can be varied depending on the signal-to-noise ratio (SNR). Second, apply a 
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1 55 dead-time correction to the photoncounting (PC) signals. For PC at high counting rates, a second 

1 56 pulse arriving at the discriminator before it has recovered from the previous pulse will not be 

1 57 counted - a period known as dead time [51]. Experiments with a function-generator-driven LED 

1 58 determine this time to be 10 ns for the high-altitude channel and 4 ns for the low-altitude channel. 

1 59 Our results show that the system dead time obeys a nonparalyzable model following a simple 

160 relationship, Eq. (1) [52], between the true count rate, C T , and measured count rates, C M , 


161 allowing the impact of dead time, Td, on the data to be removed. 


r 

C = 

T 1 -C T 

1 ^ M 1 d 


1 63 Third, remove background counts due to PMT dark counts and the sky background. Both of these 

1 64 counts are constant over all ranges and can be derived using data bins for which there are no laser 

165 signal returns. Fourth, glue the parallel analog and PC signals together [53]. The analog signal 

1 66 associated with the use of LICEL has a small time delay relative to the PC signal [54], We found 

1 67 this offset to be about 250 ns for our system by carefully comparing returns derived with clouds on 

1 68 both the analog and PC channels. The glued signal must first account for this offset. The glued 

1 69 region requires that the ratio of PC to analog signal is constant. Ratios that are not constant suggest 

170 either an incorrect background subtraction or wrong dead- time correction. The gluing threshold of 

171 the PC signal is typically 20 MHz for the Hamamatsu PMT employed in our low-altitude channel 

172 and 20-30 MHz for the EMI PMT used on the high-altitude channel. Because DIAL retrievals 

1 73 depend on the quality of both 285 and 291 signal, we glue the PC and analog signal approximately 

174 at the same altitude for both lasers to minimize the retrieval error due to gluing. Examples of the 

1 75 ratio of PC-to-analog signal and their gluing region for 285-nm signal appear in the Figure 3. The 

176 gluing threshold is 20 MHz for both altitude channels. The fifth step involves smoothing the 
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177 counts to reduce random noise. Our configuration currently employs a 5-point (5x150 = 750m) 

178 running average applied to returns from all altitudes; smoothing reduces the effective vertical 

179 resolution to 750 m. 

1 80 After initial processing, an exponential-fit correction removes signal-induced bias (SIB) 

181 from the signal returns. This bias, caused by intense light returns from near range, (also called 

1 82 signal-induced noise, SIN) appears as a slowly decaying noise source superimposed on the normal 

183 returns. The causes of the SIB are related to the regenerative effects such as dynode glow, 

184 after-pulsing, glass-charging effect, shielding effect, and helium penetration [55], SIB varies 

185 widely with different PMTs. For our case, the SIB of the EMI 9813 is larger than for the 

186 Hamamatsu 7400. SIB can persist for several hundred microseconds and can exert a strong 

1 87 influence on data at the lidar ’ s upper range where both signal and noise counts become comparable. 

1 88 With uncorrected SIB, the raw signal falls off more slowly at higher altitudes resulting in lower 

1 89 retrieved ozone values. SIB usually has more influence on the shorter wavelength channel, which 

1 90 falls off more rapidly with altitude. Unless a mechanical shutter physically blocks the optical path 

191 to the PMT to eliminate SIB, a model must characterize its behavior. Cairo et al. [56] and Zhao [57] 

1 92 have successfully used a double exponential function for this purpose. However, this correction 

193 increases measurement uncertainties because both the scaling and exponential lifetimes are 

1 94 difficult to determine without additional independent measurements. A more practical technique is 

195 to employ a single exponential fit to the residual background [39, 40, 58]. For the high-altitude 

196 channel, the function’s coefficients are automatically determined using a single exponential 

197 least-squares fit to data acquired approximately from 100 to 160 ps after data acquisition starts 

1 98 where the SIB becomes dominant. The start and length of the exponential fit varies with different 

199 channels (either wavelength channels or altitude channels), atmospheric structures, and lidar 
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200 configurations because these parameters affect the intensity of the detected signal. For our 

201 low-altitude channel, the SIB is weaker than the high-altitude channel because of the different 

202 PMT and weaker signal. However, it is difficult to automatically determine the fitting function for 

203 the low-altitude channel signal using the least-squares fitting method especially for the 285 signal 

204 because the far-range signal after background correction is not completely characterized by an 

205 exponential function (Figure 3 (b)). It is useful to optimize the exponential fitting function for the 

206 low-altitude channel using previous retrieval data and comparing the data with coincident 

207 ozonesonde profiles. The slope of the logarithm of the SIB fitting function remains for a particular 

208 configuration (i.e., outgoing power) and could slightly change for different configurations. Those 

209 retrievals corrected using the empirically-derived exponential function agree with ozonesonde 

210 profiles up to 5 km within 5% bias. Figure 3 shows the typical effect of the SIB correction and the 

21 1 comparison of the fully corrected signal and model for 285-nm signal. The model simulation 

212 employs the coincident ozonesonde measurement assuming no aerosol. 

213 3.2 DIAL retrieval 


214 The differential technique reduces the ozone retrieval uncertainty resulting from 

215 interfering aerosols and non-ozone gases. Excellent discussions concerning the DIAL technique 

216 occur in the publications by Measures [59], Kovalev and Eichinger [60], and Browell [36]. The 

217 average ozone number density, n (r+Drl2) , between range r and r+ Dr can be expressed as the 

218 summation of the signal term, n s (r+Dr/2) , the differential backscattering term, D n b (r+Drl2) , and the 

219 differential extinction term, D n e (r+Dr/2) , as following equation: 


V+Dr/2) 


>+Dr/2) 


T>n e 

Lyn (r+T)rl2) * 


221 One can write the discrete forms of the three terms at right side as follows: 
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where the subscript ‘on’ and ‘off represent the on-line (285 nm) and off-line (291 nm) 
wavelengths, respectively; P is the detected photon counts; b is the total backscatter coefficient; 

a is the total extinction coefficient excluding ozone, and D.v 03 is the differential ozone 
absorption cross-section. P , b , and a are dependent on r and wavelength. Strictly speaking, 
Ds G3 is a r dependent, as well, because it is a function of temperature, which varies with r . By 
ignoring the differential scattering and extinction from non-ozone species, the DIAL equation 
reduces to only n s . D n b arises from aerosol differential backscattering. D n e consists of 
differential Rayleigh extinction, aerosol extinction, and non-ozone gaseous absorption including 
O 2 , SO 2 and NO 2 . Measurements from a meteorological sounding can usually correct Rayleigh 
effects. We correct the aerosol effects when they are significant enough, especially in PBL. The 
aerosol correction discussion appears in Section 3.4. 

3.3 Joining retrievals from two adjacent altitude channels 

Final retrievals result from joining the data from two altitude channels with a weighted 
average. We choose to join the final ozone retrievals instead of the raw signals because the SNRs 
of the two altitude channels at the joining altitude are significantly different. If the retrievals 
derived from two different channels are statistically independent, the best estimate of these 
measurements is the 2-channel weighted average: [61] 
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where « is the ozone retrievals of channel i and the weights w. are the inverse squares of the 
corresponding statistical uncertainties ( s u , which will be discussed in Section 5), 

w,=\/s { ?. (7) 


The uncertainty of n best is 
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Typically, the low-altitude and high-altitude channels join between 3.3 and 4.4 km. Figure 
4 shows an example of a joined ozone profile as well as the combined 1 -sigma statistical 
uncertainties. 

3.4 Aerosol correction 

In a polluted area, aerosols can be a dominant error source in the lower troposphere. Based 
on Eq. (4) and (5), The vertical gradient of aerosol backscattering determines D n b and the 
magnitude of the differential aerosol extinction coefficient determines D« e . The largest aerosol 
correction usually occurs in an inhomogeneous aerosol layer (i. e., the top of the PBL). One can 
solve for the ozone and aerosol profiles simultaneously with only two wavelengths by assuming 
appropriate Angstrom exponents and constant lidar ratios [62, 63]. If a third wavelength is 
available, and is close to the DIAL wavelength pair, one can use the dual-DIAL technique [64, 65] 
to reduce the error due to aerosol. When the third wavelength is far from the DIAL wavelength 
pair, one can use the method suggested by Browell [36] to correct the aerosol interference. 
Without the third wavelength, we employ an iterative procedure to retrieve ozone and correct 
aerosol effects. To illustrate this method, start from the equation to solve for the ozone using only 
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the 291-nm signal [59]: 
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where cr 03 is the ozone absorption cross-section; and fif r) are molecular and aerosol 


» 


backscatter coefficients at range r , respectively; a^ +&r/2) and o' A r+Ar/2) represent the average 


'-(r+Ar/2) 


molecular and aerosol extinction coefficients respectively between range r and r+Dr . The 
subscript 291 is omitted for brevity, because all backscatter and extinction parameters correspond 
to 291 nm. Assuming the lidar ratio (aerosol extinction-to-backscatter ratio), S, is known for the 
291-nm signal: 

S = a A !p A , (10) 


Solving for [i A r) , Eq. (9) becomes: 


p r 2 (B M + B A 1 
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For S constant with range, a priori j3 A ef) at reference range, and downward numerical calculation, 


only one unknown variable, cc A r+Ar/2) , remains in Eq. (1 1). If we further assume: 


oc A « cc f 

tx (r+Ar/2) w ( r+Ar ) 
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is zero at the reference range. The (i A r) can be solved by Eq (11) starting from 


the 


reference range and calculating downward. Then, the first estimate for /3 A r) profile is substituted 


back into Eq. (11) to get the second estimate. However, we can use a more accurate form for 


a (r+Ar/2) aS> 
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where fif represents the value from the first estimate. With several iterations, we can get a stable 
solution for . The number of iterations required for a stable solution depends on the range 
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resolution of the signal. D n b {r+Dr/2) and Dri h ._ Dr/2) can be approximated [36] as: 
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where h is the Angstrom exponent, A/l is the wavelength separation, and B (r) is the 
aerosol-to-molecular backscatter ratio at the off-line wavelength defined as: 


B (r) = bf/bf . 


(16) 


Then, the ozone retrieval after the aerosol correction substituted into Eq. (11) recalculates the 
aerosol profile at the off-line wavelength. 

The lidar ratio (S) exhibits a wide range of variation with different aerosol refractive 
indices, size distributions, and humidity. The S measurements have been made most frequently at 
308 nm [66] and 355 nm [67-69], The S for our DIAL wavelengths was assumed to be 40 sr' 1 
constant over the measurement range. The Angstrom exponent (rj) is often seen as an indicator 
of aerosol particle size: values greater than 2 correspond to large particles like sea salt and dust, 
and values smaller than 2 correspond to small smoke particles. Most of the reported 77 for 
tropospheric aerosol are measured at wavelengths longer than 300 nm with a variation from 0 to 
2 [70, 71]. Considering h could be relatively small when it is applied in the UV region, we assume 
h = 0.5 at our DIAL wavelengths for urban aerosols [72], 
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Figure 5 shows the assumed aerosol, molecular, and ozone extinction profiles at 291 nm 
used in a model calculation to investigate the aerosol correction in the DIAL retrieval under an 
extremely large aerosol gradient condition. The hypothetical aerosol profile includes three basic 
regimes: homogeneous, increasing, and decreasing extinction. The aerosol extinction 
coefficients are arbitrarily set equal to zero below 1.2 km to investigate the downward 
propagating error. The molecular extinction profile derives from the 1976 U.S. Standard 
Atmosphere [73]. The hypothetic ozone contains a constant number density of 
1.5xl0 12 molec cm" 3 and an absorption cross-section of 1.24xl0“ 18 cm 2 molec" 1 at 291 nm. 

Figure 6 shows the comparison of the ozone retrieval both with and without aerosol 
correction as well as the calculated aerosol profile at 291 nm. This example calculation assumes 
that 7 = 0.5 and S=40 sr' 1 are known exactly, and there is no signal measurement error. With a 
range resolution of 150 m, two total iterations from Eq. (11) to (15) produce the final 
aerosol-corrected ozone retrieval. In the process of calculating the aerosol profile, 10 iterations 
of Eq. (1 1) and (13) produces a stable aerosol solution. The derived aerosol backscatter is slightly 
lower than the model at the lower altitudes in the homogeneous aerosol region. This bias, in turn, 
produces a slightly higher ozone retrieval in the corresponding region. The aerosol-correction 
procedure reduces the retrieval errors from ±60% to about ±5% . The residual errors are due to 
the numerical integration and the approximation of Eq. (14) and (15). The quality of this iterative 
procedure depends on the choice of S and h . According to Eq. (1 1), (14), and (15), S affects the 
aerosol profile retrieval while h affects only the final ozone correction. 

Figure 7 shows the sensitivity test for S and h in the aerosol correction assuming that S=40 
and h = 0.5 are the correct values. Inaccurate estimates of S or h can yield retrieval errors up to 
about 20%. Larger h will overestimate Dn e , which produces less ozone and vice versa, h has a 
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smaller impact on D n b relative to D n e due to the 4- h factor. The impact of S is larger in the 
inhomogeneous aerosol layer than in the homogeneous layer. The peak error is larger for 
underestimated S relative to overestimated S [74], 

We summarize the iterative procedure as follows: 

1. Calculate the first estimate of the ozone concentration from Eq. (3); 

2. Substitute the first estimated ozone into the Eq. (1 1) to derive the aerosol backscatter profile 
for the off-line wavelength and iterate to obtain a stable solution with Eq. (13); 

3. Calculate the differential aerosol backscatter and extinction corrections to obtain an second 
estimate ozone using Eq. (14) and (15); 

4. With the second ozone estimate, go back to step 2. 

4. Measurements 

Figure 8 shows an ozone DIAL retrieval for 12 consecutive hours from 12:56 local time, 
Aug. 9 to 03:45, Aug. 10, 2008 with 10-min temporal integration (12000 shots) and 750-m vertical 
range resolution using the data processing described in the previous section. The aerosol 
time-height curtain (Figure 8a) exhibits moderate aerosol activity below 2 km with expected 
diurnal PBL variation and shorter time-scale fluctuations due to PBL processes. The maximum 
aerosol correction in Figure 8b corresponds to an ozone adjustment of 3-4 ppbv and occurs 
between 1.5 and 2.5 km for the largest vertical backscatter gradient. There was a short duration 
cloud at 2 km at 14:00. The retrievals for the two altitude channels overlap between 3.3 and 4.4 km 
to produce the final ozone profiles (Figure 8c) that agree well with the co-located ozonesonde 
(EN-SCI model 2Z with unbuffered 2% cathode solution) launched at 13:49 local time. The 
time-height curtain of ozone’s evolution shows a very interesting structure of multiple ozone 
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346 layers in the lower atmosphere that varies with time. One can see the buildup and decay of various 

347 layers throughout this 12-hour period. The high-frequency variation in the high-altitude channel 

348 ( > 6 km) results partly from lower SNR and higher uncertainty of the SIB correction, both of 

349 which increase with altitude. Figure 9 shows the mean ozone profile and 1 -sigma standard 

350 deviation for all the 10-min vertical profiles in Figure 8 as well as the coincident ozonesonde 

351 measurement. The high-altitude channel has a standard deviation increasing with altitude due to 

352 the statistical error distribution. Its standard deviation is less than 13 ppbv below 8 km and 

353 increases to about 45 ppbv at 8.5 km where the 285 laser does not have sufficient SNR for ozone 

354 retrieval; therefore, we terminate the retrievals at 8 km in Figure 8. The standard deviation of the 

355 low-altitude channel retrievals is less than 5 ppbv below 4 km and reaches 8 ppbv at 5 km due to 

356 lower SNR. The standard deviation at 2 km is a little larger than the surrounding altitudes possibly 

357 because of larger ozone fluctuation related to the PBL top. The two altitude channels have 

358 consistent mean retrievals in the overlap region with discrepancies less than 5 ppbv and similar 

359 standard deviations at 3.3 km which most likely reflect the true ozone short-term variation 

360 considering the physically credible structure in Figure 8. The mean retrievals agree with the 

361 ozonesonde measurement within about 10 ppbv and have higher biases at the upper altitudes. 

362 5. Error Analysis 

363 We divide the error budget of the DIAL retrieval into four categories: 1) Statistical 

364 uncertainties, s x , arising from signal and background noise fluctuations; 2) Errors, s 2 , associated 

365 with differential backscatter and extinction of non-ozone gases (O 2 , SO 2 , NO 2 , etc.) and aerosols; 

366 3) Errors, s 3 , due to uncertainties in the ozone absorption cross-section; and 4) Errors, s 4 , related 

367 to instrumentation and electronics. s x is a random error; e 2 , e 3 , and s 4 are systematic errors. 
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368 


s x can be written as [38]: 



375 area of the telescope’s primary mirror, and Pl is the number of the emitted laser photons. A r must 

376 be chosen large enough to produce an acceptably small error. Figure 10 shows the estimated 

377 statistical errors for the high-altitude and low-altitude channels for a 1 0-min integration and 750-m 

378 range resolution. s x is typically less than 10% below 4 km for our low-altitude channel and could 


379 be 20% at 5 km. This altitude performance gives us sufficient overlap for the two altitude channels 

380 under most atmospheric conditions. In the high-altitude channel, e x exceeds 25% of the retrieval 


381 ozone near 8 ± 1 km, at which error altitude we terminate the retrieval. 

382 s 2 includes the interference from O 2 , SO 2 , NO 2 , air molecules, and aerosols. Table 2 


383 summaries the potential errors in the DIAL retrieval for 285 and 291 nm wavelengths due to 

384 non-ozone absorption gases. The calculation of the oxygen dimer (O 2 -O 2 ) interference includes 

385 some uncertainties due to the absorption cross-section measurement. The O 2 -O 2 absorption theory 

386 has not been entirely established [75]. The local SO 2 and NO 2 profiling data are not available. 

387 However, the estimated error due to either SO 2 or NO 2 using the latest ground observation is less 

388 than 1%. The impact caused by differential Rayleigh extinction results in an inaccuracy of less 
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than 1% using balloon ozonesonde retrievals of atmospheric density or by employing 
climatological models. 

The main concern comes from the aerosol interference, which depends on both the 
wavelengths and wavelength separation. Although the aerosol optical properties could be 
retrieved from a third wavelength, the differential effect for a DIAL wavelength pair still has some 
uncertainty due to the assumption for lidar ratio and Angstrom exponent. Within the PBL, where 
the statistical errors are small, differential aerosol backscattering and extinction dominate the error 
sources [36, 38, 40], However, it is reasonable to believe that the error due to aerosol interference 
is smaller than 20% after the aerosol correction as shown in section 3.4. 

The uncertainty in the Bass-Paur ozone cross-sections is believed to be less than 2% 
[75-77]. will be less than 3% after considering the temperature dependence. 

s 4 could be caused by a misalignment of the lasers with the telescope FOV, imperfect 

dead-time, or SIB correction. Dead time distorts the near-range signal, and SIB distorts the 
far-range signal. Because the dead-time behavior is reliably characterized, the error caused by SIB 
usually is larger than the dead-time error. For the 10-min integration data, s 4 is believed to be 

<5% at 1-4 km for our low-altitude channel and <10% for our high-altitude channel below 8 km 
based on the LED test results and the analysis of our previous data such as Figure 8 and 9. A 
summary of the errors in the DIAL measurements appears on Table 3 for a constant tropospheric 
ozone of 60-ppbv, 750-m vertical resolution, and 10-min integration. 

Figure 1 1 shows a comparison of 12 lidar retrievals and their single coincident ozonesonde 
measurement between 1300 and 1400 local time except for the first profile on Aug. 17, 2008 
(upper right panel), which was taken at 0800. The aerosol correction was made at altitudes 
between 1 and 3 km. Figure 12 shows the mean percentage differences and their standard errors of 
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412 the mean for all those retrievals. The lidar retrievals of the low-altitude channel agree with 

413 ozonesonde measurements within 10% from 1 to 4 km. The relatively high errors at about 2 km 

414 possibly relate to residual aerosol correction errors around PBL height. The lidar retrievals from 

415 the high-altitude channel agree with ozonesonde to within 20% below 8 km. The statistical error 

416 and the uncertainty associated with the SIB correction result in larger errors for the high-altitude 

417 channel above 6 km. 

418 6. Conclusion and Future Plans 

419 The RAPCD ozone DIAL system measures tropospheric ozone profdes during both 

420 daytime and nighttime using the 285/29 1-nm wavelength pair. The low-altitude receiving channel 

421 makes ozone measurements at altitudes between 1-5 km using a 10-cm telescope and 

422 Hamamatsu R7400U PMTs. The high-altitude channel measures ozone between 3 and about 8 km 

423 using a 40-cm telescope and EMI 9813 PMTs. Model calculations demonstrate that the iterative 

424 aerosol correction procedure significantly reduces the retrieval error arising from differential 

425 aerosol backscatter in the lower troposphere where the quality of the aerosol correction depends 

426 on the accuracy of the a priori lidar ratio and Angstrom exponent. A comparison of the lidar 

427 retrievals and coincident ozonesonde measurements suggest that retrieval accuracy ranges from 

428 better than 10% after the application of an aerosol correction below 4 km to better than 20% for 

429 altitudes below 8 km with 750-m vertical resolution and 10-min integration. Error sources include 

430 statistical uncertainty, differential scattering and absorption from non-ozone species, uncertainty 

431 in ozone absorption cross-section, and imperfection of the dead-time and SIB corrections. The 

432 uncertainty in the SIB correction and the statistical errors dominate the error sources in the free 

433 troposphere and could be reduced by increasing the integration time or reducing the range 

434 resolution. 
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435 Future improvements will overcome two major limitations of the current system by: 1) 

436 extending observations into upper troposphere by replacing the current transmitters with more 

437 powerful ones and shifting the current wavelengths to longer ones to make higher-altitude 

438 nighttime measurements; 2) Minimizing aerosol interference in the lower troposphere by adding 

439 a third wavelength (dual-DIAL technique). This lidar with expected improvements will provide a 

440 unique dataset to investigate the chemical and dynamical processes in the PBL and 

441 free-troposphere. The spatio-temporal variance estimates derived from the ozone lidar 

442 observations will also be useful for assessing the variance of tropospheric ozone captured by 

443 satellite retrievals. 
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711 Figure 1 . Transmitter diagram. 
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Figure 2. Diagram of the receivers and detectors. 
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722 



723 Figure 3. An example of gluing and SIB correction for 285-nm signal. The 10-min averaged data 


724 occured at 13:00 local time on Oct. 18, 2008. (a) The normalized ratio of photoncounting (PC) to 


725 analog after background and dead-time corrections for the low-altitude channel signal, (b) 
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726 Comparison of the non-SIB-corrected, SIB-corrected signal and model as well as the SIB fitting 

727 function for the low-altitude signal. The model uses the coincident ozonesonde measurement 

728 assuming no aerosol. The SIB fitting function ( exp(-l .3 -alt- 2 • 1 0” 1 ) ) uses previously retrieved 

729 data and coincident ozonesonde measurements, (c) Same as (a) but for the high-altitude channel. 

730 (d) Same as (b) but for the high- altitude channel. The coefficients of the SIB fitting function result 

731 from an empirical single exponential least-squares fit to the signal acquired from 100 to 160 ps 

732 after data acquisition starts. 

733 
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736 (a) (b) 

737 Figure 4. An example of a joined ozone retrieval for the lidar data in Figure 3. (a) Separate 

738 retrievals of the two altitude channels. The error bars represent the 1 -sigma statistical uncertainties. 

739 The gray envelope represents ±10% uncertainty of the coincident ozonesonde profde. (b) The 

740 joined DIAL retrieval from the two altitude channels and its combined 1-sigma statistical 

741 uncertainty. 
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Figure 5. Aerosol, molecular, and ozone extinction coefficient profiles at 291 nm for a model 
calculation of extreme aerosol effects. 
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751 Figure 6. Comparison of the simulated ozone retrieval without aerosol correction against that with 

752 aerosol correction using the iterative procedure. The Angstrom exponent (rj) and lidar ratio (S) 

753 were assumed to be exactly known at 0.5 and 40 sr' 1 respectively for the aerosol correction. The 

754 aerosol correction dramatically improves the ozone retrieval. 

755 
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03 (1 0 12 molec cm 3 ) 

Figure 7. Ozone retrieval using different Angstrom exponents ( 77 = 0 , 0.5, 1) and lidar ratios 
(S=20, 40, 60) in the aerosol correction. 
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(c) 


766 


767 Figure 8. Ozone DIAL retrievals made on 9-10, Aug. 2008. (a) Calculated aerosol extinction 

768 coefficient at 291 nm. The feature at 2 km; 1400 is a cloud, (b) Aerosol correction for ozone DIAL 

769 retrieval, (c) Ozone DIAL retrieval after aerosol correction. The retrieval was made with a 750-m 

770 vertical range resolution and a 10-min temporal resolution. The co-located ozonesonde marked by 

771 a triangle launched at 13:49 local time. 
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779 



Figure 9. Mean ozone mixing ratio and 1 -sigma standard deviation for the 10-min vertical profiles 
between 12:56 and 15:06 local time (excluding one cloudy profile at -1400) in Figure 8. The 
co-located ozonesonde launched at 13:49 local time. Large error bar (~45%) at 8.5 km identifies 
the high-altitude limit of retrievals (8 km). 
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780 Percentage error (%) 

781 Figure 10. Estimated statistical errors for the high-altitude and low-altitude channels using 10-min 

782 integration and 750-m range resolution. The nighttime and daytime statistical errors are modeled 

783 by using the annually averaged local ozonesonde profde, the 1976 U.S. Standard Atmosphere, an 

784 urban aerosol model [78] and the lidar parameters in table 1. The ozone profile used for summer 

785 daytime errors is assumed 20% higher than the annual average. 

786 
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795 Mean percentage difference (%) 

796 Figure 12. Mean percentage differences (dark) and their estimated 1-sigma standard error of the 


797 mean (gray) for the data in Figure 1 1 . 




799 


Table 1. Characteristics of the RAPCD ozone DIAL system 


System 

Specification 

Transmitter 

Pump lasers 

Nd:YAG, 20 Hz repetition rate, 5-7 ns pulse length, 300 mJ pulse' 1 at 
1064 nm, 50 mJ pulse' 1 at 532 nm 

Dye 

Rhodamine 590 and 610 

Emitted UV 

4 mJ pulse' 1 at 285 nm, divergence^ mrad 
3 mJ pulse' 1 at 291 nm, divergence<l mrad 

Tuning range 

277 to 303 nm for the final UV output 

Receiver 

High-altitude channel Low-altitude channel 


Telescope 


Band-pass filter 


Detector 


Signal processing 


Newtonian, 40-cm diameter, f/4.5, Welch Mechanical Designs 

1.5-mradFOV Cassegrain, 10-cm diameter, 

f/2.3, 4.3-mrad FOV 

Center wavelength at 286.4 nm with a 1 1-nm FWHM. Transmittance is 
35% at 285 nm and 20% at 291 nm 

Electron Tubes 9813QA, about 28% Hamamatsu R7400U-03, about 

quantum efficiency 20% quantum efficiency 

LICEL Transient Recorder (TR40-80), 250-MHz maximum 

photoncounting rate, 12-bit and 40-MHz analog-to-digital converter, 
25-ns range resolution 


800 

801 
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802 


Table 2. DIAL retrieval errors due to non-ozone absorption gases 


Gases 

A < 7 , differential 
absorption cross-section 
(cm 2 molec" 1 ) for 285 and 
291 nm 

References 
for Act 

Mixing ratio 
(ppbv) 

References for 
mixing ratio 

O 3 retrieval 
error (%) 

0 3 

1.15xl0” 18 

Bass and 
Paur 1981 
[76] 

60 



o 2 a 

4.5 xlO -27 

Fally et al. 
2000 [79] 

2 . 1 xl 0 8 


1.5% 

S0 2 

-4.8 xlO -20 

Rufus et al. 
2003 [80] 

13 b 

NREM 2006 
[81] 

-0.9% 

N0 2 

-2.25 xl0“ 2 ° 

Bogumil et 
al. 2003 [82] 

18 c 

NREM 2006 
[81] 

-0.6% 

Total 





±1.5% 


803 

804 

805 

806 


a due to O 2 -O 2 

b maximum 24-hr average in 1994. Latest local monitoring data available. 
c Annual arithmetic average in 1993. Latest local monitoring data available. 


50 



807 


808 Table 3. Summary of the errors in RAPCD ozone DIAL measurements* 


Errors 

Low-altitude 

channel High-altitude channel 


(1-4 km) 

(3-8 km) 

1 . s x , statistical error 

<10% 

<25% 

2. s 2 , interference by non-ozone species 



Aerosol 

<20% 

<5 % 

Non-ozone absorption gases 


<1.5% 

Rayleigh 

<1% using local radiosonde profde 

3. <s\ , due to uncertainty in Acr 03 


<3% 

4. s A , due to SIB and dead-time 

<5% 

<10% 

Total RMS error 

<23% 

<28% 


809 * The errors are estimated by assuming a 60 ppbv constant ozone mixing ratio in the troposphere 

810 for data with a 750-m vertical resolution and 10-min integration. 

811 
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